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ABSTRACT 



Analytic expressions for distance-redshift relations which have been corrected 
for the effects of inhomogeneities in the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) mass density are given in terms of Heun functions and are used to 
illustrate the significance of inhomogeneities on a determination of the mass 
parameter Q m and the cosmological constant A. The values of these parameters 
inferred from a given set of observations depend on the fractional amount of 
matter in inhomogeneities and can significantly differ from those obtained by 
using the standard magnitude-redshift (m-z) result for pure dust FLRW models. 
As an example a determination of Q m made by applying the homogeneous 
distance-redshift relation to SN 1997ap at z — 0.83 could be as much as 50% 
lower than its true value. 

Subject headings: cosmology: theory - large-scale structure of universe 
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1. INTRODUCTION 

When attempting to evaluate the mass parameter Q m and/or the cosmological constant 
A, observations of quantities such as magnitude, angular separation, and redshift are made 
on objects distant enough for curvature effects to be detected. As an example, for Type la 
supernovae (SNe la) corrected magnitudes and redshifts (m-z) are measured, plotted, and 
compared with theoretical m(O m , A; z) curves computed for the FLRW models flPerlmutter 



|et al. 1997| , |Perlmutter et al. 1998| , |Garnavich et al. 1998| ). In spite of the fact that 



the FLRW models contain only homogeneously and isotropically distributed perfect fluid 
gravity sources, one of these models is assumed to represent the "large scale" geometry of 
the universe. Relations like m(Q m , A; z) are also commonly assumed to be valid, on average. 
This latter assumption may well be incorrect for some distant observations including SNe 
la, but even if technically correct may not be useful in determining Q m and A. In particular 
if the underlying mass density approximately follows luminous matter (i.e., associated 
with bounded galaxies) then effects of inhomogeneities on relations like m(Q m , A; z) must 
be taken into account. The majority of currently observed SNe la are not being seen 
through foreground galaxies and whether or not this is due entirely to selection (rather than 
statistics) is not important. If the objects observed do not have the average FLRW mass 
density po m their foregrounds then the FLRW m-z relation does not apply to them (see 
[Kantowski 1998|) . Ultimately some SNe la should exist behind foreground galaxies ( [Rauch 



|1991|) and for these, m-z should be computed using the lensing formulas. These formulas 



( [Bourassa fc Kantowski 1975| and |Cooke fc Kantowski 1975| ) contain source-observer, 



deflector-observer, and source-deflector distances, respectively D s ,Dd, and Dd s , all of which 
depend on the mass density in the observing beam, excluding the deflector. These distances 
will not be given by the standard FLRW result if the observing beam contains less than the 
average FLRW mass density. 
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In §|2] the average area-redshift equation ([]]) for a light beam traveling through a FLRW 
Swiss cheese universe is given and its solution is related to the luminosity distance-redshift 
relation D^z). In §|] the solution of this equation is given for the case where gravitational 
lensing can be neglected. The new result of this paper, Dp(z) without lensing for FLRW 



Swiss cheese can be found in equations fl42| ) and ([47]), and for the special case f2 = 1 in 



equations ( |A2| ) and ( |A3| ) of Appendix A. In §^ numerous m-z plots are given to illustrate 



the effects of inhomogeneities and some conclusions are drawn. It is argued that if 
homogeneities are not taken into account when attempting to determine fl m and A, errors 
as large as 50% could be made. Even though the Di(z) given here has been derived using 
the exact Swiss cheese cosmologies, the result are valid for observations in essentially any 
perturbed pressure-free FLRW models in which lensing can be neglected. 

Inhomogeneous models of the Swiss cheese type and their associated optical equations 
discussed here are often mistakenly attributed to Dyer and Roeder (see Appendix B). 

Appendix C contains some useful simplifications for evaluating the real-valued Heun 
functions needed in the analytic m-z relations given here. Appendix C also contains six 
useful lines of Mathematica code which numerically evaluates and plot these same m-z 
relations. 



SWISS CHEESE OPTICS 



Some years ago the author ( [Kantowski 1969| ), used the "Swiss cheese" cosmologies to 



study the effects of local inhomogeneities in the FRW mass density on the propagation 
of light through an otherwise homogeneous and isotropic universe. That analytic work 
was undertaken because prior results computed in perturbed FRW models were suspect. 
In particular, results were inconsistent with weak lensing results, e.g., on average, the 
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luminosity of a distant object was not given by the FRW result ( [Bertotti 1966| ). Numerous 
sources of error were suggested but particularly mistakes inherent in using perturbative 
gravity were suggested. For example the FRW relations between radius, redshift, and 
affine parameters were (and still are for approximate GR solutions) assumed valid in the 
presence of mass perturbations. At that time |Kantowski 1969| put to rest any question 



of the possible existence of an effect on the mean luminosity; theoretically it could exist! 
Because the Swiss cheese models are exact solutions to the Einstein equations the accuracy 
of the FRW relations between radius, redshift, and affine parameters could be directly 
established. Today Swiss cheese itself is under attack as the source of the 'erroneous' 
prediction (see [Frieman 1997| and [Wambsganss et al. 1997| and additionally see the 



related work of Premadi et al. 1998| ). The intent of this particular paper is not to defend 
Swiss cheese predictions against these attacks, that can be done elsewhere. The reader 
should not dismiss distance-redshift predictions made by Swiss cheese models because of 
their non-physical distribution of matter, e.g., cheese and holes. Because these models 
contain the only two types of gravitational curvature (Ricci and Weyl) that affect optical 
observations, and because they are fairly flexible in including density perturbations, they 
should adequately describe optical observations at least as far as z = 1 (a distance to which 
galaxies and other inhomogeneities are thought to have undergone only minor changes). 
The conjectured extension of the validity of the optical equation (P used here beyond the 
Swiss cheese models has been made frequently since it was first derived and is argued by 



Schneider et al. 1992 in their Sections 4.5.2 and 4.5.3. 



The purpose of this paper is to extend analytic results for distance-redshift relations 
in inhomogeneous FRW models to FLRW, i.e., to include the cosmological constant (see 
Kantowski et al. 1995| and |Seitz fc Schneider 1994j ) . [Kantowski 1969| derived a 2nd order 



intergal-differential equation [see (43a) of that paper or (1) of |Kantowski et al. 1995| as well 



as the equivalent 3rd order differential equation, (43b) of that paper or (6) of [Kantowski et 
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al. 1995] for the average cross-sectional area A of a beam of light starting from a distant 
source and propagating through a A = Swiss cheese universe (see Figure 1). The solution 
of this equation, with appropriate boundary conditions, gives all average quantities relating 
to distance-redshift. This equation and its derivation were easily extended to include a 
cosmological constant by Pyer fc Roeder 1974| . However, as will be seen below, extension 



of the equation's analytic solutions of |Kantowski et al. 1995| and [Seitz fc Schneider 1994 



is somewhat involved and the special functions required are much less familiar to the 
math/physics community. Weinberg sign conventions will be again used (|Misner et al. 
19731) . 



As a light beam from a distant SN la propagates through the universe (Figure 1) the 
cheese of the model produces the same focusing effect as does the transparent material 
actually appearing within the beam. The holes in the cheese with their condensed central 
masses reproduce the optical effects of the remaining Friedmann matter that has been 
condensed into clumps, e.g., galaxies. The extended equation for the average area A 
traversing the universe, that is randomly focused by numerous clouds of transparent matter 
and lensed by numerous clumps is: 

y/A A* 2 p iZm [i + Z) ' [i) 
where prime (') is differentiation with respect to an affine parameter, 

d 



' = -(1 + zfsjl + Q m z + n A [(l + z)-* - 1] — , (2) 
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{f) **B.IU.[*£±#d,. (3) 

A On J Z 



and (£ 2 ) is the average of (a/A) 2 , the square of the wavefront's shear over its area,Q 

15 Pi u o f A 2 (l + z) 6 
Po 

In (|l|) pd (D is for dust) is the average mass density of all transparent material interior to 
light beams used to observe the given objects and pi (I is for inhomogeneous) is the average 
mass density of all types of clumpy material systematically or statistically excluded from 
the light beams. The average shear term in (|I|) comes from the Weyl (conformal) curvature 
tensor of inhomogeneous material exterior to the beams. The Pd/po term comes from the 
Ricci tensor of transparent material within the beams. For the SNe la observations, pr, 
would certainly include those ubiquitous low mass neutrinos (if they exist) as well as other 
transparent material not confined to galaxies, while pi would contain all matter clumped 
with galaxies. If there is no correlation of mass and light, deciding what goes in p^ and 
what goes in pi is problematic, and the relative value becomes another unknown parameter 
of the theory. The current Friedmann mass density is the total po — Pd + Pi and the 
curvature parameter Qn = Q m + Q\ consists of a mass part and a cosmological constant 
part: 



n ^ = ^7TT2- and n ^ = ^m- ( 4 ) 



8nGp Ac 

mf and n ^3H 2 

Inclusion of the cosmological constant A (using FLRW rather than FRW) only modifies the 
functional relationship between redshift and affine parameter (|^), see Pyer fc Roeder 1974 . 



The unitless gravitational lensing parameter B is defined in equation (A2) of |Kantowski 



et al. 19951 and its effects on the solution of equation ([!]) are described in pyer fc Roeder 



1974 . In this paper analytic solutions to equation ([!]) will be given for Bq = 0, i.e., for 



2 The form of (£ 2 ) depends on structure details of the clumps. What is given in (3) is for 
objects completely condensed into opaque spherical masses. This particular type of Swiss 
cheese clumping is expected to produce maximum lensing. In the following sections we will 
be interested in observations where even maximum lensing is negligible. 
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distance-redshift in any Swiss cheese model where (£ 2 ) is neglected. In [Kantowski et al. 



1995 it was argued that even maximal lensing effects (B ^ 0) are not expected to be 
significant when observing SNe la at z < 1 and, as pointed out above, when lensing events 
do occur, the observed magnitudes should be analyzed by using the lensing formula and not 
by incorporation into the m-z relation. The resulting B = equation (|II]) represents the 
equation for the average area of a light beam only Ricci focused by part of the mass density, 
Po)- Such a light beam is not conformally lensed by inhomogeneities (pi = po — po) 
that remain exterior to the beam. If Weyl lensing is infrequent, a distribution of areas 
will occur for which the Bq = equation gives the maximum value for the area (i.e., the 
lower bound on the distribution of luminosities). In |Kantowski 1998| the resulting m-z 



is appropriately dubbed the 'intergalactic' magnitude-redshift relation because it is m-z 
without galactic focusing. If significant galactic lensing is an unusual event as apparently is 
the case with SNe la beams passing exterior to galaxies, the 'intergalactic' m-z approximates 
the 'mode' value (the most likely) .0 If galaxies are compact (20 kpc) the 'intergalactic' m-z 
relation should be more useful in determining Q m and fl\ from SNe la observations than 
is the mean m-z relation (standard FLRW relation). If galaxies are more diffuse (200 kpc) 
exact modeling of the lensing galaxies will be important. 

To relate the differential equation (|I|) to observations, consider a source at redshift z s 
radiating power SV into solid angle SQ. The flux received by an observer at z = in area 
A\ is given by T = 5V/A\ (1 + z s ) 2 . The two factors of (1 + z s ) can be thought of as 
coming separately from the redshift of the observed photons and their decreased rates of 



3 



This assertion is consistent with that part of the numerical work of [Holz fc Wald 



|1998| which treated galaxies as condensed objects. When galaxies were modeled by 200 
kpc isothermal spheres the 'mode' moved towards the 'mean' and away from the minimum 
(intergalactic) value as expected of a more homogeneous model. 
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reception. The definition of luminosity distance is motivated by this result, i.e., 



(5) 



The observed area A\ is evaluated by integrating equation (|l|) from the source z = z s to 
the observer z = with initial data which makes the wave front satisfy Euclidean geometry 
when leaving the source (area=radius 2 x solid angle): 



<A\ 

<u/a\ 



0. 



(6) 



dz ""^H s (l + z s y 

where in FLRW the value of the Hubble parameter at z s is related to the current value H 
at z = by 



h s = h (i + z s ) y/i + n mZs + n A [(i + Zs )- 2 - 1]. 

The series solution of equation (0) , combined with (||) and @ is: 



(7) 




1 + f2 A — 



1 

+ 2 



+ 



2^m f ^ ~ ~~ ~~ ' ~~ ' " ' 



10 - 2v(v + 1) - 5a 



-B t 



' 6 



Ha 5 + - 



5 + z/(z/ + l) 



Q m + — fi^ + 5fi A - -Oa 



4 + 3ft r 



+ Ofz 5 



where the source redshift z s has been simplified to z and pj / po has been replaced for later 
convenience by a clumping parameter v, < v < 2, 



1 + 24(p//po) - 1 = v(y + 1) 

Po 



2 fh; 6 

This series is useful for understanding the low-redshift sensitivity of Di to the various 
parameters; e.g., Q m and Oa appear in the z 2 term, the clumping parameter v first 



(9) 



(8) 
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appears in the z 3 term, whereas the lensing parameter B doesn't appear until the z 4 term. 
Additionally, analytic results computed in the next section can be checked by comparison 
with this series. 



3. THE ANALYTIC SOLUTION FOR D e (z) WHEN LENSING CAN BE 

NEGLECTED 

In this section the general B = solution of ([!]) will be given for boundary conditions 
appropriate for De(z). If apparent-size (angular) distances are desired the reader has only 
to compute D < (z) = Di(z)/(\ + z) 2 . The new solution appears in (|42|), (|47D , (|A2|), and 
( |A3| ) expressed in terms of Heun functions HI. All previously known special solutions are 
limiting cases of the general solution ( f47l) [see (fL8D and Appendix B]. To solve equation ([[]) 
it is first rewritten as: 



;i + zfyji + n mZ + q a [(i + z )- 2 - 1] X 

d x o i ~ — — : — 7 d 



;i + zfyjl + Q m z + Q A [(1 + z)- 2 - 1] T y/A(z) 



az v dz 

4 

This equation is often attributed to Dyer-Roeder ( Pyer fc Roeder 1972| , Pyer fc Roeder 



f (3 + i/)(2 v) ^ + g)Sy fi£j = Q (1Q) 



1973 ) in the literature, (see Appendix B for some history of this equation). To date 



only numerical solutions have been obtained when 7^ 0, e.g., see |Asada 1996| , |Suto 



[Matsubara 1996| , [Kayser et al. 199"7| . It can be put into a recognizable form by changing the 



independent variable from z to y and the dependent variable from J A(z) to h 

y = y (l + z)= ^ m (1+z), 




h = 0- + Z )\i rrV (11) 
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The resulting equation is 

d 2 h | + dh i y („ + i) y + i h _ Q 

dy 2 y 3 + y 2 - b n dy y 3 + y 2 - b n 
where b n = -VL 2 m tt K /(l - fi m - fi A ) 3 . When the cubic y 3 + y 2 ~b n = (y - y x )(y - y 2 ){y - y 3 ) 
is factored, (|i~2|) simplifies to a recognizable form of the Heun equation (see |Ronveaux 1995| , 
Erdelyi 1955| , |Whittaker fc Watson l927HHeun 1889| ): 



d 2 h ( 7 | 8 | e \dh ^ a(3 y - q h — (13) 

d?/ 2 V-V2 V-Vz) dy (y - yi){y - y 2 ){y - y 3 ) 

where (|12|) requires 

7 = 5 = 6= 2' 
1 

a = 

^ = + 1), 

q = 1, (14) 

and additionally the three roots to be constrained by: 

2/12/22/3 = b n = -Vt 2 m Vt A /(l -Vt m - Q A ) 3 , 
2/i + 2/2 + 2/3 = -1, 

2/i2/2 + 2/i2/3 + 2/22/3 = 0. (15) 



The Heun equation is slightly more complicated than the hypergeometric equation; it 
possesses four regular singular points in the entire complex plane rather than just three. 
In the form given by (|H|) one of the two exponents of each finite singular point (y±, y 2 , 2/3) 
vanishes and the other exponent is given respectively by (1 — 7, 1 — 5, 1 — e).[] The point at 

4 Recall that an exponent gives the analytic behavior of a solution within the 
neighborhood of a regular singular point, e.g., h = (y — y 1 ) 1_7 (l + ci(y — yi) + ■ ■ •) ■ 
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oo is the fourth singular point and its exponents are a and (3. For the point at oo to also 
be regular (i.e., for this to be a Heun equation) all exponents must sum to a value of 2, 
equivalent ly 

a + /5 + l = 7 + 5 + e. (16) 

For fll^D this necessary constraint is satisfied. From flip] ) it follows that at least one root 
has to be real and complex roots must come in conjugate pairs. For convenience y\ will be 



chosen as real throughout. This Heun equation (|T3|) is conveniently expressed in terms of a 
Riemann P-symbol as: 



P { 



Vl 2/2 2/3 oo 
a y q 
1-7 1-5 1-e (3 



> . 



(17) 



The first 4 columns of ( \17\ ) are the 4 regular singular points and their 2 exponents, and 
the 5th column is the independent variable, all analogous to the Riemann P-symbol for 
the hypergeometric equation. The 6th column is the constant q from the numerator of 
the coefficient of h in the Heun equation [when put into the standard form of fljjP]. The 
hypergeometric equation is uniquely specified by information about its regular singular 
points but Heun requires the extra parameter q. Because the three finite singular points 



of this Heun equation have values of 1/2 for their nonvanishing exponents, (jl2|) can be 
transformed into the Lame' equation. In this paper, solutions of (|T2|) will be given as 
local Heun functions and in a following paper they will be expressed as Lame' functions. 
When A = (=>- bn = 0) equation ( |I~2"D has only three regular singular points i.e., this 
Heun equation simplifies to an equation of the hypergeometric type. Additionally the 
corresponding Lame' equation reduces to the associated Legendre equation. Solutions for 
A = can be written either as combinations of hypergeometric functions or as associated 
Legendre functions [see (|i~8D and (|B7|) below]. To motivate the form of the £7a ^ solution, 



13 



the fl A — solution will be given first [[Kantowski et al. 1995| and [Seitz fc Schneider 1994 



Dg.(VL m = Q , Vl A = 0, u; z) 
c 1 

x 



„ ,,,„/, (V V 3 3 ^ \ / V V 1 1 1 — Qn 

1 + tt z ) + 1 2F1 - + 2, - + -; v + -; 1 - Q ) 2 Fi( 1, \--v\ 

v U; \2 '2 2' 2' / \ 2 '2 2'2 ' 1 + Q z 



[1 + zf 



;i + fi ^) 3/2+ " /2 



v v 1 1 \ /z/ z/3 3 1 - Q \ 

1, ; - - v. 1 - fl ) 2 F 1 I - + 2, - + -; u+ -; 

2 ' 2 2'2 ' V V2 '2 2' 2'1 + lW 



(18) 



The expected form of the Q A 7^ solution follows the above where the hypergeometric 
functions 2 F\ are replaced by local Heun functions HI [see (|42|) and 



From (|T|) the three singular points (2/1, 2/2, 2/3) are chosen from the six permutations of 
the three roots: 



F 3 



^+ 
1 

1 + — + u_ 

f_ 

— i- 
e 3 is 
1 + + e l 3<u + 



(19) 



where 



-1 + b + Jb(b - 2 



with 



l-b + Jb(b-2 



b=^-b n = -^tf m tt A /(l -tt m - tt A ) 



(20) 



(21) 



The locations of the three finite singular points (2/1,2/2,2/3) are determined by the value 
of the single parameter b, (—00 < b < 00). Some important values are shown as contours 
in Figure 2. 
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The standard form for the Heun equation ordinarily has its singularities at (0, 1, a, oo). 
The simple linear transformation 



c 



y-yi 
2/2 — 2/1 ! 



moves 



2/i 

2/2 
2/3 



0. 
1. 



2/3-2/1 2/i( 2 + 32/i) (Vs-Vif 



2/2-2/1 (2/2 — 2/i) 2 2/i(2 + 3yi) 
The latter two equalities are consequences of the useful identity: 



(2/2 - 2/i) (2/3 - 2/i) = 2/i(2 + 3y 



1), 



which results from (|15"1). In terms of the new variable (, ( |i~2"|) becomes 



d 2 /i 1/11 1 



and the value of q changes to: 



The new Riemann P-symbol is: 



C(C-l)(C-a) 
l + M^ + l)2/i 



2/2-2/1 



(22) 



(23) 



(24) 



(25) 



(26) 



Pi 



1 a 00 

-2 c 



2/2 -l/l 



111 £+1 

2 2 2 2 



> . 



(27) 



See the figures in Figure 3 for locations of a in the complex plane and the trajectories of 
C(z) starting with ( (the value of ( at zero redshift [see ([□]) and (0)]) for the following 
three cases: 

b < — ► y 1 = real, y 2 = 2/3, and |a| = 1, 
< b < 2 — >■ 2/1 , 2/2, 2/3, a are all real, 

2 < b — ► yi = real, yi = f/3, and \a\ = 1. (28) 
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Figure 3 gives the proper choices for the three roots (2/1, 2/2,2/3) from the six possible 
orderings of (Yi, Y 2 , Y 3 ) in each of the three b domains. It also contains values for a, q, the 
new variables (, and Co- Hyperbolic and trigonometric variables, £ and <p, can be used to 
parameterize the values of the three roots (rather than b) and they are also given in Figure 
3. 

Boundary conditions on h come directly from its definition (|TT|) and the desired 
boundary conditions on \fA [see 



'A 
dy/A 



KCs) = 0, 



dz 



HJ1 + z s 



dh 

dz 



c 



Equation (|5|) then relates Di to the value of h at the observer, 

D E (z s ) = (l + z,)h(z = 0). 
Using these boundary conditions on two independent solutions hi & h 2 of 



hi(( s )h 2 (( s ) - hiQhiCs) 



m 


IdC 






1 dz 


Z s - 



where h = ^. From (|D and (|TT|) 



l/o(l + z)-y 1 

V2-V1 

yo-yi 



Co 

V2 - yi 
d( yo 



dz y 2 - yi ' 

The denominator of h(() in ([!]]) can be evaluated using the Wronskian of 



(29) 



(30) 



gives 



(31) 



(32) 



C(C-l)(C-a) 



(33) 
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where Cw is a constant. The square root in this term can be evaluated using 



\/c(c-D(c-<.) = v f TW ^[ s 



3/2 



y 3 + y 2 -b n = (i + z)^=^i + n m z + n A [(i + z)-2-i}. (34) 



With from fl30|) and # s from (|7|), equations ([31]) and (|33|) give the desired result: 

= - 77 + n MOfcCCo) - fc(C.)MG>)] • (35) 

#oV"m(j/2 - yi) 1/2 (CV) 

Figure 2 shows domains in the (fi m , Oa) plane separated by 6 = 2, 6 = oo, and \(o\ = 1. 

The b = oo contour is equivalent to Oo = 1. These contours are important because they 

separate domains for which different choices of the two independent solutions hi(() and 

/2.2(C) must be taken. The \( \ = 1 contour divides the \( \ < 1 domain where solutions 

about the singular point £ = are chosen from the \(q\ > 1 domain where solutions about 

oo are chosen. These choices are necessary for convergence of the local Heun functions. For 

the Qa = case, only analytic expressions about oo were required [see (fHf )]. 

The Heun Function Hl(a, q; a, /3, 7, 5; () is the analytic solution of fl39|) defined by the 
infinite series ([36|), see [Ronveaux 1995| . It converges in a circle centered on ( = which 
extends to the nearest singular point 1 or a. This solution is analogous to the 2 F\ solution 
of the hypergeometric equation but unfortunately does not appear in any of the common 
computer libraries. When cq = 1 is chosen (as will be done here) the series is: 



00 



Hl(a,q;a,P,j,5;() = 1 + ^cvC, (36) 

r=l 

where the c r are constrained by a three term recursion relation (take c_i = 0): 

P r c r _i - (Q r + q)c r + R r c r+1 = 0, (37) 

with 

P r = (r-i + a )(r-l + p), 
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Q r = r[(r - 1 + 7)(1 + a) + a6 + e] 
Rr = (r + l)(r + 7)a. 



(38) 



The e parameter is not included as an argument in ii7(a, g; a, /3, 7, 5; () because of the 
constraint (0). This series corresponds to the zero exponent for the regular singular point 



( = and will be taken as hx(C) m dH) when |£ Q | < 1. The second independent solution is 

h 2 (C) = C^Hl(a, q H - a n , p n , 7//, 5; {), (39) 



where four parameters have changed, 



Qn = (a8 + e)(l -7) +q, 
ct n = a + 1 — 7, 
P„ = /? + 1 - 7, 
7// = 2- 7. 



(40) 



The constant CV in the Wronskian can be evaluated for the ( ~ expansion using /ii 
and above, 

CV = (41) 
This gives an expression for the luminosity distance appropriate for \( Q \ < 1, 

c(l + z) 



D e (tt m ,tt A ,v;z) 



2/o(2/o ~ 2/i) 

i?o|V^^i(2 + 3 yi ) 



x 



HI a. 



1 + \u{u + l) yi 
2/i(2 + 3yi) 



z/ z/ + l 1 1 |/o (l+z) -yi 
yi(2 + 3yi) 



2' 2 '2'2' 



3 + (z/ 2 + z/-3)yi ^ z/-l z/ + 2 3 1 yo-2/i 
xiii a, = — va; 



4^2/1(2 + 33/! 



2 ' 2 '2'2' 



2/1(2 + 37/1 



ly (l + z) -yi„, I 3 + (u 2 + u-3) yi v-l is + 2 3 1 y (l + z) - Vl 

-til a, ===== — va; z — , — z — , -, -; — — va 



xi?Z a, 



y°-yi v 4^(2 + 3^ 



2/1(2 + 32/0 



2 2 '2'2' ^(2 + 3^ 
v v+1 1 1 2/0-2/1 



a ' 2' 2 '2' 2' 



2/i(2 + 32/i 



(42) 
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where the source redshift z s has again been replaced by z. The required values of y\ and 
a can be found in Figure 3 for all three b domains and yo = Q m /(1 — Qq) is the value of y 
at z = [see ( PD ]. Even though the above Heun functions contain complex arguments and 
parameters, they are real valued functions of the real redshift varable z. As soon as these 
functions become available in Mathematica, expressions fl42|) and (f47|) will be immediately 
useful. Untill then, simpler expansions suitable for \a\ = 1 are indicated in Appendix C. 

The solution similar to ([42]) but suitable for \( Q \ > 1 is given by choosing 




(43) 



and 




(44) 



where seven parameters have now changed: 




(3 = -(3 + 5 + e 



7 = 1 + a - (3 




a + 1 - 7, 




/3 + 1-7, 



In 



2- 7 . 



(45) 



For this choice of hi and \i2 the constant in the Wronskian (|33|) becomes 



[3 + a = - 




(46) 



From (|35|) the resulting expression for the luminosity distance is 
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(S2 m , S2 A , i>; z) = x 

#o(z/ + §)>AW i +z)- yi/yo 



u+l 

y (l + z) - yA 2 
2/o — 2/i / 

1 4 - u 2 + i/(l - 2v)y x 1 i/l-i/ l-2i/l Jyjp + ^Vi) 1 



°' 4,/yi(2 + 3?/i) v 7 "' 2 ' 2 2 ' 2' y (l + z) - y x y/a 



xHI\- (3 + ^)(l-i')-(3 + 2i/)(l + % 1 1 1 + J/ z/ + 2 3 + 2z/ 1 _ 0/i(2 + 3y a ) l 



40/i(2 + 3yi) v/a' 2 2 2 '2' y - 2/i 

v_ 



w (i + «) -yi 



,„.l (3 + i/)(1-i/)-(3 + 2i/)(1 + %i 1 l + z/z/ + 2 3 + 2z/l Jyi(2 + 3 yi ) i 
x/it — 



4^1(2 + 32/1) v/a' 2 2 ' 2 ' 2' ?/o(l + ^) -2/1 ^ 

(47) 



1 4 - z/ 2 + i/(l - 2i/)j/i 1 1/ l-v l-2v l_\/yi(2 + 3yi) l 



v a ' 4^1(2 + 32/1) v 7 ^' 2 ' 2 2 '2' 2/0-2/1 v 7 ^ 

The special case of f2o = fi m + — 1 can be obtained from fl42|) and (|47|) by taking the 
appropriate limits. Some details of this process along with the resulting luminosity distance 
are given in Appendix A. For those values of Vt m and Q\ where \(q\ < 1 it is clear that for 
large enough values of z, |£| > 1 and hence (f42"D is no longer valid (HI no longer converges). 
For some values in the (fi m , Oa) plane above the \Co\ — 1 contour in Figure 2, fl42"l) will not 
converge for a SNe la range of z ~ 0.5, but for most values it does. 

In the next section several plots of magnitude vs. redshift are made to illustrate the 
importantce of take clumping into account when attempting to determine Q m and Q\. 
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4. m-z PLOTS FOR CLUMPY UNIVERSES & CONCLUSIONS 

In this section several magnitude-redshift plots are given to illustrate the effects that 
density clumps can have on the m-z relation and consequently on a determination of Q m and 
A made by using this relation. Because m-z depends differently on Q m and as a function 
of redshift for the FLRW models, both parameters could in principle be determined from 
a sufficient quantity of accurate SNe la data. Clumping provides an additional parameter 
v which complicates any such determination. As can be seen from (|8|) the dependence of 
m-z on this additional parameter could also be determined by enough data. However, such 
a triple determination is certainly more complicated. What will be done here to illustrate 
the effects of the v parameter is to plot multiple m-z curves for various values of all three 
parameters u, Q m , and Q\. In all plots the unit of distance is taken to be c/H Q . In these 
figures De is plotted on a magnitude scale, 5 Log De (i.e., the distance modulus plus 5 Log 
10pcHo/c). In Figure 4, fl\ is held fixed while v and fl m are varied and in Figure 5, Q m is 
held fixed while v and Q\ are varied. In Figure 6, = Q m + £7a = 1 is fixed while all three 
parameters vary. 

In Figure 7 the sensitivity of observed magnitudes to variations of fl m is illustrated by 
fixing z = 0.83 and Q\ = 0.1. In Figure 8 a similar plot is given showing the sensitivity 
to variations of fl^- The importance of the clumping parameter is easily seen from these 
last two figures. If the distance modulus of a source such as SN 1997ap at z = 0.83 were 
precisely known (e.g., see the two sample horizontal lines in Figure 7) then a determination 
of Q m could be made, assuming Q\ were somehow known. Likewise, from Figure 8, a 
determination of Q\ could be made if fl m were somehow known. From Figure 7 the reader 
can easily see that the determined value of fl m depends on the clumping parameter v. The 
Q m value will be about 95% larger for a v = 2 completely clumpy universe than it will be 
for a v = completely smooth FLRW universe. Equivalently, fl m could be underestimated 
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by as much as 50% if the FLRW is used. The maximum underestimate is reduced to 33% 
at the smaller redshift of z = 0.5 (see a similar result for Q\ = in [Kantowski et al. 1995|) . 
These conclusions are not sensitive to the value of Q\. 

Slightly different conclusions follow from Figure 8 about The discrepancy in the 
determined value of Q\ is A£7a ~ —0.14 for v — 2 compared to v — 0, and is not sensitive 
to the distance modulus. The discrepancy is halved, A^a ~ —0.07, at a smaller redshift of 
z = 0.5. 

A minimal estimate of the quantity of data required to begin distinguishing between 
the various v values can easily be made. At z — 0.5 the differences in observed magnitudes 
of a SN la in a v = (100% smooth FLRW) and a v = 2 (100% clumpy) universe is about 
Am ~ 0.02 if Q m ~ 0.2, and Am ~ 0.09 if Q m ~ 0.8. These differences are not sensitive to 
A. With corrected-intrinsic and observed magnitude uncertainties of ±0.2, |Branch 1998|, 



data on over 200 SNe la will be required if we live in a low density universe and over a 
dozen if we live in a higher density one. 

The results presented here (fil2"|),(|4"?|),( |A2| ), and ( |A3| ) for the 'intergalactic' distance- 
redshift relation are quite general. They contain corrections (for mass inhomogeneities) 
to the standard FLRW result, applicable to observations where gravitational lensing can 
be neglected, i.e., observations where the conformal (Weyl) curvature doesn't produce 
significant average shear in (JTJ) . Even though the original area equation (P was rigorously 
established for a particular type of Swiss cheese model, the resulting equation which neglects 
lensing fllPf ) is expected to be widely applicable to observations at redshifts of z — 1 and 
less. Application of its solution to a given set of observations requires that the average 
fraction of the mass density contained in the observing beams (i.e., the v parameter) be 
determined. This fraction obviously depends on the number as well as the type of object 
observed. Collecting CMB radiation at wide angles is likely to produce a v = value but 
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observing a few dozen SNe la might well result in a value close to v = 2 (i.e., we might in 
fact live in a universe where mass, dark or otherwise, is primarily associated with galaxies). 
If a significant fraction of the universe's mass density is clumped on galactic scales, then 
the effects of these clumps on SNe la observations should be taken into account by using 
the lensing formulas rather than by decreasing v to zero. Recent numerical work by [Horz fc 



Wald 199"E| confirms the assertion that, given galaxy clumping, the cross-sectional area of 
a typical light beam will not follow the FLRW area-redshift relation. Instead the area will 
follow more closely one of the 'intergalactic' m-z relations given here, until a lensing event 
occurs. The new luminosity distances presented here represent the theoretical minimum 
of the observed magnitudes and are especially applicable to situations where lensing is 
infrequent (i.e., where the most probable value is closer to the min than the mean). Because 
[Holz fc Wald 1998| did not include any diffuse transparent matter, the applicable m-z 



relations given here are those with v = 2. For A = and Qq < 1 it is the Dyer-Roeder 
solution ( |B3|) and for A = and Qq — 1 its the v = 2 solution of Dashevskii & Slysh (|B^ ). 



The v = (standard FLRW) result represents the theoretical 'mean' for m-z for a 
universe in which only weak-lensing events occur. For extremely non-symmetric probability 
distributions, the "mean" is not likely the best estimator - in this case the "most probable" 
is likely better, [Sivia 199T) . 
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A. THE SPECIAL CASE: fi = VL m + Q f 



A complete derivation of the Qq = Q m + — 1 case can be done by introducing a new 
independent variable in fll2|), y — > — y/A where A = f2 m + — 1 and then taking the limit 
A — ► 0. The resulting differential equation which replaces (|12|) has the same exponents 
given in ( |14|) but has q = 0. The three finite regular singular points are now located at 
(2/1, 2/2, 2/3) — (^m^A) 1 ^ 3 ( — 1) e~ i7r//3 , e J?r / 3 ). When the modified equation fll2|) is transformed 
by (P^D a modified (EH) results which is described by the Riemann P-symbol: 



P< 



1 e i7r / 3 00 



000-1 



11 1 v+X 
222 2 



V3 



4^3 



> . 



(Al) 



The solutions to this simpler Heun equation with boundary conditions appropriate for 
luminosity distance Dg can be obtained directly from ( [12] ) and ( fj7|) by simply substituting 
?/o = ^m/A, yi = — (fi^A^^/A, a = e^ 3 , and then taking the limit A — > 0. This gives 
an expression for the luminosity distance, appropriate for |£ | < 1, 



A?(fi m ,fi A ,^; z) 
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The expression for the luminosity distance appropriate for \( \ > 1 is obtained by applying 
the limiting proceedure to fl47|), 
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These expressions were used to produce Figure 6 of The Q A — > limit of (|A3|) results in 
the solution (|B2l) below, first given by [Dashevskii fc Slysh 1966 . 



B. PREVIOUSLY KNOWN SOLUTIONS FOR D e (z) WHEN B = 

Until now, analytic solutions to the average area equation (P, neglecting lensing (i.e., 
putting (£ 2 )=0), have been found only for A = 0. The earliest solutions were written down 
before the equation was formulated by [Kantowski 196S . The standard homogeneous FRW 



solution was given by [Mattig 1958| . It is the pi = (i.e., v = 0) solution of ([Tj|) and 

Dt{n m = n Q , n A = o, v = o- z ) = |n ^ + (n - 2) (yji + n z - 1) } . (bi) 
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The Q Q = 1 solution was given by Dashevskii fc Slysh 1966| . [] 



Hn (u + ±) 



2 ' 2 



(B2) 



The = 2 (i.e., po = 0) solution is due to Pyer fc Roeder 1972| : 



D e (Q m ,Q A = 0,u = 2;z) 
c QJl + z) 2 
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This result can be rewritten using the identity 



1 + z) 
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sinh 1 



no(i + 2) ~ 2 ^v^-v/r^, 



(B4) 



in the form actually given by pyer fc Roeder 1972| . When fio > 1 equation 



is analytically continued using y/1 — Qq — > ±iy/Qo ~ 1, which simplifies by using, 
sinh -1 (is) = isin _1 (x) to give a form containing only real variables, 
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(B5) 



Zcl'dovich 1964| had given the Q m = 1, A = 0, v — 2 solution along with the first 



derivation of (|i0|) restricted to that particular case. Zel'dovich seems to be the first to 
recognize the importance of inhomogeneities on FRW distances and to attempt to make the 



needed modifications, pashevskii fc Zel'dovich 1965| later extended the equation to fl m 7^ 1 
but only for v = 2. It was pashevskii fc Slysh 1966| that produced the first v 7^ 2 or 0, 
A = 0, Q m 7^ 1 version of (|T0| ) along with the fl m = 1 solution ( |B2| ). 
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The above solution was the first analytic solution to the A = version of ([!]) and cited 
Kantowski 1969| as the origin of the equation. This solution represents the most extreme 
effects that clumps perhaps have on the standard Mattig result (|B1|) and is therefore one of 
the more interesting and useful solutions. For this solution no matter converges the light 
beam. The following solution ( [B6D appeared in [Dyer fc Roeder 1973| and after the above 
solution but did not cite [Kantowski 1969| . Instead, the derivation was essentially repeated 
neglecting lensing (i.e., putting the shear term to zero) and assuming, without justification, 
that the FRW relation between redshift and affine parameter was valid. Many authors have 
subsequently dubbed ([!]) with B = as the Dyer-Roeder equation. As pointed out in the 
previous footnote several other non-rigorous derivations of FRW versions of (pi]) already 
existed by 1966. Besides making frequent and successful use of the equation, Dyer- Roeder 's 
contribution to its development was to give two special solutions to it and to extend the 
equation to include A ^ 0. It is the authors opinion that if anyone deserves to have their 
name attached to a curved-space optics equation it's R. K. Sachs because equations such 
as ([!]) are direct applications of [Sachs 1961| . The v = 1 (i.e., pn = ^Pi) solution of ([!]) that 
appeared in Pyer fc Roeder 1973| was: 



D e (Q m ,Q A = 0, v= l;z) 
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The general A = solution ( |T8[ ) was only recently obtained and can be written using 
associated Legendre functions ([Kantowski et al. 1995| and [Seitz fc Schneider 1994|) as 
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The general A ^ solutions fl%2]) , P?|) , ( |A2| ),and ( |A3| ) take on a form similar to this except 
the associated Legendre functions P 2 & Ql are replaced by Lame' functions. 
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C. A SIMPLIFICATION FOR HI AND MATHEMATICA CODE FOR 

D e (tt m ,tt A ,v; z) 

Because the Heun functions are not yet available as standard computer routines, some 
programming skills are required if expressions ([42| ), (fl7]) , ([A"2]) , and (|A3|) are to be made 



used of. For those interested in computing the Heun functions used in these expressions the 
modified form of fl36|) used to produce Fig's. 4-8 is given and for those simply interested in 
plots of the luminosity distance, a few lines of Mathematica code are also given. 

For most points in the Q m -Q\ plane of Figure 2, HI as given by ([36]) contains complex 
coefficients C T as well as a complex variable (. Because the needed HI are in fact real 
functions of a real variable z it is convenient to modify this series to make its coefficients 
real. The modification is slightly different for expansions about Q = and ( = oo. For the 
expansions about ( = a convenient coefficient to iterate is C r defined by: 

CrC = u s 92rf , ° r , ( -Ay - yiY, (Ci) 

and for the expansion about oo, 

Cr (}) = m& (^) ■ (C2) 

The new C r are real polynomials and hence easily evaluated and simplified. 

For those who want only to plot m-z relations a Mathematica routine which 
numerically integrates ( |12"D with boundary conditions @ appropriate for obtaining 
(H /c) x Dg(Q m , Qa, v\ z) follows. As input the program requires values for Om, 01, and nu 
which stand for Q m , Q\, and v respectively. 
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yO= 0m/(l-0m-01) 
bO= -Om"2*01/(l-Om-01)"3 
P= (l+3/2y)y/(y"3+y"2-bO) 
Q= -(l+nu*(nu+l)*y/4)/(y"3+y"2-bO) 
distance= NDSolve[{d" [y]+P*d'[y]+Q*d[y] == 0, 
d[y0]==0, d'[y0]==l/y0}, d, {y yO, 2*y0}] 
Plot[Evaluate[(l+z)*d[yO*(l+z)] /. distance], {z, 0, 1}] 



This routine was used to check for errors in (M2J) , fl47|) , (|A2|) , and 
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Fig. 1. — Radiation beam of cross- sectional area A propagating through a Swiss cheese 
universe from distant source to observer. 

Fig. 2. — The fl m -fl\ plane and various domains required for the needed Heun functions. 
Fig. 3. — This table/figure contains various parameters needed for evaluating the luminosity 



distance ^(fi m , £Ia, v\ z) as given by ([42]) and (f47|). All parameters are functions of the 
cosmic parameters Q m , and the clumping parameter v. Three columns are given for 
the three domains of the parameter b separated by contours in Figure 2. The figure in 
each respective column gives the locations of the regular singular point a in the complex 
£-plane, the orbit of ( in (|32|) as a function of redshift z, including the starting point ( . 
Domains in the Q m -Q\ plane that correspond to |£o| < 1 and |£o| > 1 are shown in the small 
reproductions of Figure 2 included in each column. 

Fig. 4. — Magnitude-redshift relation, 5 log 10 ^Di(VL m , Vt^ = 0.1, u; z), as a function of 
redshift z for three values of Q m and three values of v. 

Fig. 5. — Magnitude-redshift relation, 5 log 10 ^fD([Vl m = 0.2, ft\, v; z), as a function of 
redshift z for three values of fl\ and three values of v. 

Fig. 6. — Magnitude-redshift relation, 51og 10 ^Di(Q m , Q\ = 1 — Q m , v\ z), as a function of 
redshift z for three values of Q m and Qa (^o = + = 1) and three values of v. 

Fig. 7. — Magnitude, 5 log 10 ^Di{VL m , fl\ — 0.1, v; z — 0.83), as a function of fl m for three 
values of v. 

Fig. 8. — Magnitude, 5 log 10 ^fDi{VL m = 0.2, Q^, v\ z = 0.83), as a function of for three 
values of v. 
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useful parameters for finding singular points 
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